Coupled Cluster Treatment Of An Interpolating Triangle/Kagome Antiferromagnet 
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The coupled cluster method (CCM) is applied to a spin-half model at zero temperature which 
interpolates between a triangular lattice antiferromagnet (TAF) and a Kagome lattice antiferromag- 
net (KAF). The strength of the bonds which connect Kagome lattice sites is J, and the strength 
of the bonds which link the non-Kagome lattice sites to the Kagome lattice sites on an underlying 
triangular lattice is J'. Our results are found to be highly converged, and our best estimate for 
the ground-state energy per spin for the spin-half KAF (J' = 0) of —0.4252 constitutes one of the 
most accurate results yet found for this model. The amount of classical ordering on the Kagome 
lattice sites is also considered, and it is seen that this parameter goes to zero for values of J' very 
close the KAF point. Further evidence is also presented for CCM critical points which reinforce the 
conjecture that there is a phase near to the KAF point which is much different to that near to the 
TAF point ( J = J'). 

PACS numbers: 75.10.Jm, 75.50Ee, 03.65.Ca 



Our knowledge of the zero-temperature properties of 
lattice quantum spin systems has been enhanced by the 
existence of exact solutions, mostly for s = 1/2 one- 
dimensional systems, and by approximate calculations 
for higher quantum spin number and higher spatial di- 
mensionality. Of particular note have been the density 
matrix renormalisation group (DMRG) calculations j| 
for one-dimensional (ID) and quasi-lD spin systems, al- 
though the DMRG has, as yet, not been so conclusively 
applied to systems of higher spatial dimensionality. Sim- 
ilarly, quantum Monte Carlo (QMC) calculations |j| 
at zero temperature are limited by the existence of the 
infamous sign problem, which in turn is often a conse- 
quence of frustration for lattice quantum spin systems. 
We note that for non-frustrated systems one can often 
determine a "sign rule" [Q, [j| which completely circum- 
vents the minus-sign problem. 

A good example of a spin system for which, as yet, no 
sign rule has been proven is the spin-half triangular lat- 
tice Heisenberg antiferromagnet (TAF). The fixed-node 
quantum Monte Carlo (FNQMC) method || has however 
been applied to this system with some success, although 
the results constitute only a variational upper bound for 
the energy. Other approximate methods (?|, ||, ^, [l(| have 
also been successfully applied to the spin- half TAF, and 
most, but not all, such treatments predict that about 50% 
of the classical Neel-like ordering on the three equivalent 
sublattices remains in the quantum case. In particular, 
series expansion results 0j give a value for the ground- 
state energy of E g /N=— 0.551, although the correspond- 
ing value for the amount of remaining classical order of 
about 20% is almost certainly too low. This spin-half 
TAF model therefore constitutes a very challenging prob- 
lem for such approximate methods. However, the spin- 
half Kagome lattice Heisenberg antiferromagnet (KAF) 
poses an even more difficult problem, because, like the 
TAF, not only is it highly frustrated and no exactly prov- 
able "sign rule" exists, but also the classical ground state 



is infini tely degenerate. Careful finite-sized calculations 
|pd| , [l2[ [f| have however been performed for the quan- 
tum spin-half KAF, and these results indicate that none 
of the classical Neel-like ordering seen in the TAF re- 
mains for the quantum KAF model. The best estimate 
for the ground-state energy of the KAF via finite-sized 
calculations |l3| stands at E g /N=—0A3. Furthermore, 
series expansion results j?]] indicate that the ground-state 
of the KAF is disordered. Indeed, a variational calcula- 
tion which utilised a dimerised basis also found that 
the ground state of the KAF is some sort of spin liquid. 

In this article we wish to apply the coupled cluster 
method (CCM) to a model which interpolates between 
the spin-half TAF and spin-half KAF models, henceforth 
termed the J- J' model (illustrated in Fig. |l|). The 
Hamiltonian is given by 



H = J ^2 S i 



{i,k} 



s k 



(1) 



where runs over all nearest-neighbour (n.n.) bonds 
on the Kagome lattice, and {i, k} runs over all n.n. bonds 
which connect the Kagome lattice sites to those other 
sites on an underlying triangular lattice. Note that each 
bond is counted once and once only. We explicitly set 
J = 1 throughout this paper, and we note that at J' = 1 
we thus have the TAF and at J' — we have the KAF. 

We now briefly describe the general CCM formalism, 
although for further details the interested reader is re- 
ferred to Refs. |l(], The exact ket and bra ground- 
state energy eigenvectors, |^) and of a general many- 
body system described by a Hamiltonian H , 



(2) 



are parametrised within the single-reference CCM as fol- 
lows: 



|*) = e s |$) 



1*0 



2 




(a) 



(b) 



FIG. 1: The J-J' model is illustrated in diagram (a), where 
the bonds of strength J between Kagome lattice sites are in- 
dicated by the thick solid lines and the non-Kagome bonds of 
strength J' on the underlying triangular lattice sites are indi- 
cated by the "broken" lines. The triangular lattice Heisenberg 
antiferromagnet (TAF) is illustrated in diagram (b), and it is 
noted that the two models are equivalent when J = J'. The 
quadrilateral unit cells for both cases are also illustrated. The 
J-J' model contains four sites per unit cell, whereas the TAF 
has only one site per unit cell. 



(tf I = ($\Sc 



S = 1 
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(3) 



The single model or reference state |$) is required to have 
the property of being a cyclic vector with respect to two 
well-defined Abelian subalgebras of multi-configurational 
creation operators {Cf} and their Hermitian-adjoint de- 
struction counterparts {CJ = {Cty}. Thus, |$) plays 
the role of a vacuum state with respect to a suitable 
set of (mutually commuting) many-body creation oper- 
ators {Cj}. Note that Cj\$) = 0, V I ^ 0, and that 
Cq~ = 1, the identity operator. These operators are fur- 
thermore complete in the many-body Hilbert (or Fock) 
space. Also, the correlation operator S is decomposed en- 
tirely in terms of these creation operators {Ct }, which, 
when acting on the model state ({C/|$)}), create exci- 
tations from it. We note that although the manifest Her- 
miticity, = |*)/(*|*)), is lost, the normalisation 

conditions (^l^) = ($|*) = ($|$) = 1 are explicitly im- 
posed. The correlation coefficients {Si, Si} are regarded 
as being independent variables, and the full set {Si, Si} 
thus provides a complete description of the ground state. 
For instance, an arbitrary operator A will have a ground- 
state expectation value given as, 

A = = ($|S*e- s Ae s |$> =A({S I ,S I }\ . (4) 

We note that the exponentiated form of the ground- 
state CCM parametrisation of Eq. ([|) ensures the correct 
counting of the independent and excited correlated many- 
body clusters with respect to |$) which are present in 
the exact ground state |$). It also ensures the exact 
incorporation of the Goldstone linked-cluster theorem, 
which itself guarantees the size-extensivity of all relevant 
extensive physical quantities. 



The determination of the correlation coefficients 
{Si, Si} is achieved by taking appropriate projections 
onto the ground-state Schrodinger equations of Eq. (||). 
Equivalently, they may be determined variationally by 
requiring the ground-state energy expectation functional 
H ({Si, Si}), defined as in Eq. (||), to be stationary with 
respect to variations in each of the (independent) vari- 
ables of the full set. We thereby easily derive the follow- 
ing coupled set of equations, 

SH/SSi = =S> ($|C7e~ s iJe 5 |$) = 0, 7^0 ; (5) 
SH/SSi = (®\Se- s [H, C*+]e s |$) =0, I + (6) 

Equation (||) also shows that the ground-state energy at 
the stationary point has the simple form 



E g =E g ({S I }) = ($\e- s He s \$) 



(7) 



It is important to realize that this (bi-)variational formu- 
lation does not lead to an upper bound for E g when the 
summations for S and S in Eq. (|3|) are truncated, due to 
the lack of exact Hermiticity when such approximations 
are made. However, one can prove that the important 
Hellmann-Feynman theorem is preserved in all such ap- 
proximations. 

In the case of spin-lattice problems of the type con- 
sidered here, the operators Cf become products of spin- 
raising operators s^T over a set of sites {k}, with respect to 
a model state |$) in which all spins points "downward" in 
some suitably chosen local spin axes. The CCM formal- 
ism is exact in the limit of inclusion of all possible such 
multi-spin cluster correlations for S and S, although in 
any real application this is usually impossible to achieve. 
It is therefore necessary to utilise various approximation 
schemes within S and S. The three most commonly em- 
ployed schemes previously utilised have been: (1) the 
SUBn scheme, in which all correlations involving only 
n or fewer spins are retained, but no further restriction 
is made concerning their spatial separation on the lat- 
tice; (2) the SUBn-m sub-approximation, in which all 
SUBn correlations spanning a range of no more than m 
adjacent lattice sites are retained; and (3) the localised 
LSUBm scheme, in which all multi-spin correlations over 
all distinct locales on the lattice defined by m or fewer 
contiguous sites are retained. 

For the interpolating J-J' model described by Eq. (Q) , 
we choose a model state | $) in which the lattice is divided 
into three sublattices, denoted {A,B,C}. The spins on 
sublattice A are oriented along the negative z-axis, and 
spins on sublattices B and C are oriented at +120° and 
— 120°, respectively, with respect to the spins on sublat- 
tice A. Our local axes are chosen by rotating about the 
y-axis the spin axes on sublattices B and C by —120° 
and +120° respectively, and by leaving the spin axes on 
sublattice A unchanged. Under these canonical transfor- 
mations, 
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TABLE I: CCM results for the ground-state energy per spin 
and sublattice magnetisation of the TAF and KAF models us- 
ing the LSUBm approximation with m — {2, 3, 4, 5, 6}. CCM 
critical values, J' c , of the J-J' model (with J = 1), which are 
themselves indicators of a phase transition point in the true 
system, are also given. Comparison is made in the last row 
with the results of other calculations. 





KAF 


TAF 


J-J' 


m 


EJN K 


M K 


E g /N 


M K 


J' c 


2 


-0.37796 


0.8065 


-0.50290 


0.8578 




3 


-0.39470 


0.7338 


-0.51911 


0.8045 


-0.683 


4 


-0.40871 


0.6415 


-0.53427 


0.7273 


-0.217 


5 


-0.41392 


0.5860 


-0.53869 


0.6958 


-0.244 


6 


-0.41767 


0.5504 


-0.54290 


0.6561 


-0.088 


oo 


-0.4252 


0.366 


-0.5505 


0.516 


0.0±0.1 


c.f. 


-0.43°. 


0.0[a/. 


-0.551 6 . 


0.5 C . 





a See Refs. M 
b See Ref. ft 
c See Refs.lH 
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FIG. 2: CCM results for the ground-state energy per spin of 
the J-J' model (with J = 1) using the LSUBm approxima- 
tion with m = {2,3,4,5,6}. The boxes indicate the CCM 
critical points, J' , and a simple extrapolation in the limit 
m — » oo implies that J' c = 0.0 ± 0.1. 
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The model state | <3?) now appears mathematically to con- 
sist purely of spins pointing downwards along the z-axis, 
and the Hamiltonian (for J = 1) is given in terms of these 
rotated local spin axes as, 



H 



= E + ^ + s * s i - 4^ - 

J' Yl {~\ s i s * + ^( s i4 + s l s k - 44 - *i 4) 



{i^k} 

A 



+ ^(4 s k + s i 4) - y^44 + s i s k )} 



(9) 



Note that i and j run only over the Nk sites on the 
Kagome lattice, whereas k runs over those non-Kagome 
sites on the (underlying) triangular lattice. N indicates 
the total number of triangular-lattice sites, and each 
bond is counted once and once only. The symbol — > in- 
dicates an explicit bond directionality in the Hamiltonian 
given by Eq. (|J), namely, the three directed nearest- 
neighbour bonds included in Eq. (||) point from sublat- 
tice sites A to B, B to C, and C to A for both types 
of bond. We now perform high-order LSUBm calcula- 
tions for this model via a computational procedure for 
the Hamiltonian of Eq. The interested reader is 

referred to Ref. for a full account of how such high- 
order CCM techniques are applied to lattice quantum 
spin systems. 



We note that for the CCM treatment of the J-J' model 
presented here the unit cell contains four lattice sites (see 
Fig. 0). By contrast, previous calculations [[To] for the 
TAF used a unit cell containing only a single site per 
unit cell. Hence, the J-J' model has many more "funda- 
mental" configurations than the TAF model at equivalent 
levels of approximation. However, we find that those con- 
figurations which are not equivalent for the J-J' model 
but are equivalent for the TAF have CCM correlation co- 
efficients {Si, Si} which become equal at the TAF point, 
J' = 1. Hence, the CCM naturally and without bias re- 
flects the extra amount of symmetry of the J-J' model at 
this one particular point. This is an excellent indicator of 
the validity of the CCM treatment of this model. The re- 
sults for the J-J' model at J' — 1 thus also exactly agree 
with those of a CCM previous treatment of the TAF. Our 
approach is now to "track" this solution for decreasing 
values of J' until we reach a critical value of J' c at which 
the solution to the CCM equations breaks down. This is 
associated with a phase transition in the real system p0| , 
and results for J' c for this model are presented in Table~|. 
A simple "heuristic" extrapolation of these results gives 
a value of J' c — 0.0 ± 0.1 for the position of this phase 
transition point. This result indicates that the classical 
three-sublattice Neel-like order, of which about 50% re- 
mains for the TAF, completely disappears at a point very 
near to the KAF point ( J' = 0). 

The results for the ground-state energy are shown in 
Fig. H and in Table |[ These results are seen to be highly 
converged with respect to each other over the whole of the 
region < J' < 1. A simple heuristic extrapolation may 
be attempted for these results for varying J' by plotting 
LSUBm results for m = {3,4,5,6} against 1/m 2 and 
performing a linear extrapolation of this data as was done 
previously Jl(| for the TAF only. These results are given 
in Table | for the KAF and TAF models. We believe that 
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FIG. 3: CCM results for the sublattice magnetisation of the 
J-J' model (with J = 1) using the LSUBm approximation 
with m — {2, 3, 4, 5, 6}. Note again that a simple extrapola- 
tion of LSUBm critical points in the limit m — > oo implies 
that J' = 0.0 ± 0.1. 



the extrapolated results are among the most accurate 
results for the ground-state energies of the TAF and KAF 
ever found. 

We now wish to consider how much of the original clas- 
sical ordering of the model state remains for the quantum 
system. Previous calculations for the TAF |l(| took the 
average value of sf (after rotation of the local spin axes) 
where i runs over all lattice sites. However, if one does 
this one also includes non-Kagome lattice sites, and when 
J' = the spins on this site would be effectively "frozen" 
to the original direction of the model state at these sites. 



Hence, we believe that the correct order parameter for 
this model is the average value of (again after rota- 
tion of the local spin axes) where k runs only over the 
Nk Kagome lattice sites. We may thus write this as, 



M 



K 



2 



N K 

fc=i 



(10) 



The results for M K are presented in Fig. || and in Table 
H Again, we extrapolate these results for the KAF by 
plotting LSUBm results for m = {3,4,5,6} against 1/m 
and performing a linear extrapolation of this data, as 
was done previously flfiJ] for the TAF. The extrapolated 
result for the KAF point probably lies too high. However, 
the LSUB6 result goes to zero very close to the KAF 
point, and so CCM results are fully consistent with the 
hypothesis that, unlike the TAF, the ground state of the 
KAF does not contain any Neel ordering. 

It has been shown in this article that the CCM may be 
used to provide highly accurate results for the ground- 
state energy of the J-J' model (with J = 1) which 
interpolates between the TAF and KAF models. In- 
deed, the extrapolated results for the ground-state en- 
ergy for the KAF of E g /N K =-0A252 and for the TAF 
of E g /N—~ 0.5505 are among the most accurate yet de- 
termined for these models. Furthermore, the amount of 
classical ordering (evaluated on the Kagome lattice sites 
only) yields results which are fully consistent with the 
hypothesis that the KAF is fully disordered. CCM criti- 
cal points also reinforce the conjecture that the classically 
ordered phase evident for the TAF breaks down very near 
to the KAF point. 



[1] S.R. White and R. Noack, Phys. Rev. Lett. 68, 3487 
(1992); S.R.White, Phys. Rev. Lett. 69, 2863 (1992); 
S.R.White, Phys. Rev. 5 48 10345 (1993). 

[2] K. J. Runge, Phys. Rev. B 45, 12292 (1992); ibid. 45, 
7229 (1992). 

[3] A.W. Sandvik, Phys. Rev. B 56, 11678 (1997). 

[4] R.F. Bishop, D.J.J. Farnell, and J.B. Parkinson, Phys. 
Rev. B 61, 6775 (2000). 

[5] R.F. Bishop and D.J.J. Farnell, in Advances in Quantum 
Many-Body Theory, Vol. 3, edited by R.F. Bishop, K.A. 
Gernoth, N.R. Walet, and Y. Xian, World Scientific, Sin- 
gapore (2000) - in press. 

[6] M. Boninsegni, Phys. Rev. B 52, 5304 (1995). 

[7] R.R.P. Singh and D.A. Huse, Phys. Rev. Lett. 68, 1766 
(1992). 

[8] B. Bernu, P. Lecheminant, C. Lhuillier, and L. Pierre, 
Phys. Scripta, T49, 192 (1993); Phys. Rev. B 50, 10048 
(1994). 

[9] T. Jolicoeur and J.C. LeGuillou, Phys. Rev. B 40, 2727 



(1989). 

[10] C. Zeng, D.J.J. Farnell, and R.F. Bishop, J. Stat. Phys., 
90, 327 (1998). 

[11] B. Bernu, P. Lechimant, and C. Lhuillier, Physica 

Scripta. T49. 192 (1993). 
[12] P. Lecheminant, B. Bernu, C. Lhuillier, L. Pierre, and P. 

Sindzingre, Phys. Rev. 5 56, 2521 (1997). 
[13] C. Waldtmann, H.-U. Everts, B. Bernu, P. Sindzingre, 

C. Lhuillier, P. Lecheminant, and L. Pierre, Eur. Phys. 

J. 5 2, 501 (1998). 
[14] C. Zeng and V. Elser, Phys. Rev. B 51, 8318 (1995). 
[15] R.F. Bishop, Theor. Chim. Acta 80, 95 (1991). 



